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way by combining the SFR in the local universe, the re- 
sults of ground-based spectroscopic redshift surveys out 
to z ~ 1 (Lilly et al. 1996) and the results of deep photo- 
metric surveys like the Hubble Deep Field (HDF) in which 
galaxies photometric redshift have been determined by the 
Lyman dropout selection technique (Madau et al. 1996, 
hereafter M96). The resulting cosmic star formation his- 
tory has been summarized in the so-called "Madau plot" 
which shows a rapid increase up to z ~ 1 and then a 
decrease at higher z with a peak at z ~1.5 - 2. How- 
ever, one of the large uncertainties in deriving the CSFR 
from the comoving luminosity density is dust obscura- 
tion. Indeed, the faint objects have been selected optically 
and therefore may represent only the tip of the iceberg 
because dust absorbs large amounts of UV-optical rest- 
frame light, causing underestimation of the SFR by a large 
factor. This factor is still a matter of debate, for exemple 
Pettini et al. (1998) have estimated the extinction on the 
Lyman Dropout galaxies (at z > 2) and found a factor 
of ~ 3 while Meurer et al. (1997) derived for the same 
objects a factor as large as 15. Finally, there might exist 
a population of star-forming galaxies at high-z which are 
completely enshrouded in dust and therefore remain un- 
detected in the optical/near infrared as revealed by the 
ISOCAM-HDF survey (Aussel et al. 1998, Elbaz et al. 
1999 ). Moreover, new estimates of the the global star 
formation rate have been reported. The first one is de- 
rived from the analysis of deep ISO images of the CFRS 
by Flores et al. (1999, hereafter F99) who revised previ- 
ous star formation density by a factor of ~ 3 upward at 
< z < 1 and the second one comes from an analysis of 
Key words: cosmology: diffuse radiation-galaxies: evolutionLyman Break Galaxies in the HDF by Steidel et al. (1999, 



Abstract. We study the evolution of stellar content and 
the chemical enrichment of the universe averaged over 
the whole population of galaxies by means of a series 
of chemo-spectrophotometric models that take into ac- 
count the metallicity and dust obscuration effects. We 
investigate various classes of cosmic star formation rates 
(CSFR) histories consistent with current estimations. We 
are able to reproduce a variety of observational constraints 
such as the emissivities at rest-frame 0.44, 25, 60 and 100 
/im of the local universe and also the overall shape of 
the cxtragalactic background light from UV/NIR galaxy 
counts and the cosmic infrared background (CIB) from 
DIRBE/FIRAS measurements. We find that the CIB at 
140 /im is crucial for discriminating between the CSFR 
histories. The best-fit model to this constraint seems to 
favour mid-infrared derived CSFR at low-z and a flat 
CSFR at higher z consistent with the most recent estima- 
tions, although the shape of the CSFR at high redshifts 
has little impact on the FIR/submm part of extragalactic 
background. We suggest that the bulk of the CIB energy 
is produced by a population of moderately obscured nor- 
mal galaxies lying at < z < 1. We then derive the global 
chemical enrichment and stellar content of the universe 
and find that this model predicts metallicities in good 
agreement with the metallicity in DLAs given that some 
outflow of metal-enriched gas from galaxies is assumed, 
but it overproduces the current present-day stellar mass 
density and NIR luminosity density. 



ISM-intergalactic medium 



1. Introduction 



The history of the cosmic star formation rate (CSFR) 
in the universe has been inferred in an unprecedented 
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hereafter S99) who found a higher SFR at z ~ 4 than pre- 
vious estimate by Madau et al. 1998 (hereafter M98) and 
therefore found no evidence for a significant decline of the 
CSFR at high redshifts. 

The history of CSFR can also be observed by its im- 
print on the background radiation. The first detection in 
DIRBE data of the IR/submm background by Puget and 
his collaborators (1996) at a level comparable and even 
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higher than the optical counterpart, suggests that a large 
fraction of the energy of young stars is absorbed by dust 
and reradiated in the IR/submm. Understanding the na- 
ture and redshift of the sources responsible for the infrared 
background is of great importance for galaxy evolution 
studies. Several deep surveys have been undertaken in the 
submm wavelengths range at 450 and 850 /im with the 
SCUBA Bolometer Array and ISO satellite (at 175 /an), 
in order to resolve this background (Kawara et al. 1997, 
Barger et al. 1998, Hughes et al. 1998, Eales et al. 1999). 
Recent advances in high resolution spectroscopy of high 
redshift systems like damped Lya systems (DLAs), has 
provided us with precious information on early phases of 
chemical enrichment (Pettini et al. 1997, hereafter P97) 
and neutral gas content of the universe (Storrie-Lombardi 
et al. 1996) casting some light on the cosmic chemical evo- 
lution. The aim of this paper is to explore the predictions 
that models of stellar evolution self-consistently coupled 
with chemical evolution code would make on the metal and 
stellar enrichment history of the universe. Briefly, given 
the cosmic star formation history we calculate the unat- 
tenuated spectral luminosity density at different rcdshifts. 
Adopting a model of dust absorption we then compute the 
attenuated comoving luminosity density at different red- 
shifts and wavelengths which is then compared to available 
observations. Combining this with the extragalactic back- 
ground will allow us to constrain the models. Predictions 
on the cosmic chemical enrichment and on the stellar con- 
tent of the universe are then derived. 
This paper is organized as follow: section 2 describes the 
models and the inputs we used. In section 3, we discuss 
our models predictions in the light of available photomet- 
ric constraints, such as the comoving luminosity densities 
from UV/vis to FIR and the cosmic infrared background 
radiation (CIB). In section 4 we derive the chemical en- 
richment history of the universe and its stellar content. 
Throughout the paper we use a value of the Hubble con- 
stant H = 50 km s _1 Mpc -1 , a density parameter f2 = 1 
and Qa = 0. 

2. The models 

In this work we study the global properties of galaxies 
averaged in large comoving volumes. We are in particu- 
lar interested in the stellar and heavy element content of 
the universe that our models would predict using avail- 
able observational constraints. This is based on chemo- 
spectrophotomctric models that follow in a self-consistent 
way the evolution of stars, dust and metals. These mod- 
els are based on stellar population synthesis code which 
is coupled to simple chemical evolution models and dust 
absorption. The basic ingredients of our models are de- 
scribed bellow (more details can be found in Devriendt et 
al. 2000, hereafter DGS) 



2.1. The stellar emission 

Our models are based on up-dated libraries of the- 
oretical stellar tracks and spectra allowing for the build- 
ing of new sets of spectrophotometric evolutionary mod- 
els that take into account the effect of metallicity. These 
models use the Geneva group stellar evolutionary tracks 
(Schaller et al. 1992, Schaerer et al. 1993a, b) for five metal- 
licities Z = 0.001, 0.004, 0.008, 0.02, 0.04. The main stel- 
lar evolution phases include all the important stages ex- 
cept for low mass stars (less than 1.7 M Q ) those tracks 
stop at the Giant Branch tip. We thus use the more re- 
cent track evolution of low mass stars which go from the 
Zero-main sequence up to the end of EAGB for only two 
metallicities Z = 0.02 and Z = 0.001 (Charbonnel et 
al. 1996). For the other metallicities we simply interpo- 
late between the available Geneva tracks. Choosing the 
Geneva group tracks allows the use of the Z-dependcnt 
stellar yields of Maeder (1992) which insures the con- 
sistency when coupling spectrophotometric and chemical 
evolutions. For the stellar spectra we use theoretical spec- 
tra by Kurucz (1992), with the advantage that it avoids 
the use of the transformations of bolomctric luminosity log 
Lboi to magnitude My , of log T e / / to Spectral Type, and 
to luminosity classes, which are rather uncertain for the 
hottest ones. For the coldest stars (K and M-type stars) 
with T < 3750 K we prefer to use models from Bessell 
et al. (1989, 1991a,b) for M Giants (for various Z) and 
from Brett (1995a, b) for M dwarfs, with a single metal- 
licity Z=0.02. The initial helium content of the gas from 
which stars form is computed with an initial He fraction 
Y = 0.24 and assuming AY/AZ = 3. 

2.2. Modelling the Cosmic Chemical Evolution 

A model of chemical evolution allows the abundance 
evolution in the ISM and the abundance distribution in 
stars by solving a simple set of equations (Tinsley 1980). 
For a population of galaxies, the chemical evolution equa- 
tions can be expressed in terms of comoving densities of 
stars, neutral gas and heavy-elements Vt g , f2*, Viz in unit 
of the critical mass density p c = 7 x 10 10 M Mpc~ 3 : 

^L = -ZVl^ + Vl z -Vl z out (1) 

Z is the mean abundance of the ISM of galaxies defined by 
the ratio Vlz/Vl g . Vl z ut is the comoving rate of outflowing 
metal-enriched gas out of galaxies. In the case of closed- 
box models this term vanishes. 

We use the Z-dependent yields of M92 with moderate 
mass loss for massive stars. For intermediate stars we use 
the Renzini and Voli (1981) yields. We do not make the 
assumption of instantaneous recycling but we track metals 
ejected by stars at each timestep and let new stars form 
out of the enriched gas, assuming instantaneous mixing of 
the metals. 
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2.3. Dust spectra 

A non-negligible part of the energy released by stars 
in the UV/optical wavelengths is absorbed by dust and 
re-emitted in the IR/submm. A normal spiral such as 
our Galaxy emits ~ 30% of its light in the far-infrared 
(see Dwek et al. 1998 and reference therein). Taking dust 
obscuration into account is crucial for the census of stel- 
lar mass, metals and energy. Models of SEDs which cor- 
rectly connect UV/vis and IR/submm windows are there- 
fore needed. However modelling dust extinction is not an 
easy task since we have to address critical issues such as 
the dust properties, its geometrical distribution with re- 
spect to stars and its chemical evolution. Models of SEDs 
that take into account the extinction and emission effects 
of dust in a consistent way have been developed by sev- 
eral authors (Franceschini et al. 1991, 1994), Fall, Chariot 
and Pei (1996), Silva et al. (1998), Pci et al. (1999, here- 
after PFH), but all these models assume a simple relation 
between dust and metal abundance of the gas and do not 
have a detailed modelling of the starlight re-processing by 
dust grains. More detailed and self-consistent models of 
SEDs called STARDUST have recently been developped 
by DGS. Here we will follow DGS to compute the global 
SED from UV/vis to IR/submm wavelengths. 



2.3.1. Dust absorption 

The first step is to compute the optical depth. We will 
first assume a simple prescription given by: 



rx(t) 



(-7—) z q t v 
A v 



(2) 



where (A\/Ay)o is the extinction curve in our Galaxy (so- 
lar metallicity) from Mathis et al. (1983). In principle one 
should relate T\(t) to the mean H-column density which 
is proportional to the mean gas density times the gaseous 
disk scale (Nii(t) = M g {t)/iTr 2 oc p g asr g ). To account for 
this, we introduce a "fudge" factor ry. 
As in Guiderdoni and Rocca-Volmerange (1987, hereafter 
GRV87), we assume a face-on "slab" geometry distribu- 
tion where the gas and the stars which contribute mainly 
to dust heating are distributed with equal height scales in 
the disks. Following GRV87, if T\(t) is the optical thick- 
ness of the disks at wavelength A and time t, the mean 
internal extinction correction (averaged over inclination 
angle i) is then given by: 



A x (t) = -2.5 log < 



1. — exp{—a\T\{t) / cosi) 



a\T\{t)/cosi 

The factor a\ = (1 — ujx) 1 ^ 2 takes into account the 
effect of the albedo lo\. This "slab" geometry is an inter- 
mediate case between the "screen" geometry, where the 
dust layer lies in front of the stars layer and the "sand- 
wich" geometry, where the dust layer is trapped between 



two star layers, geometries which leads respectively to a 
larger and a smaller absorption (Franceschini and An- 
dreani 1995, Andreani and Franceschini 1996). 
We have also to take into account the evolution of ry 
with rcdshift. As already mentionned, Ty is in principle 
related to the mean H-column density Nh and depends 
on the geometrical structure of the gas. From scaling ar- 
gument ry oc Pgasrg we would expect that Ty cx (1 + z) a 
with a = 2. However, other factors might evolve, such 
as the properties of the dust or the actual geometry of 
the gas distribution within the disk so that the redshift- 
dependence might deviate from the simple scaling relation. 
In this case, the slope index a is taken as an adjustable pa- 
rameter whose value is determined when consistency with 
the HDF rest-frame UV luminosities is obtained. 
The above modelling of the mean optical thickness is very 
crude and probably does not reflect the complexity of dust 
extinction. In the following, we will use the GRV87 pre- 
scription in which the mean optical thickness is explic- 
itly related to the mean metallicity and to Nh- Following 
GRV87 and Franceschini et al. (1991) we will assume that 
the mean face-on optical depth of the average gas disk at 
time t and wavelength A is given by : 

MO-OS) fW( .„WL, ). <«> 



2.1 10 21 at cm" 2 



where the mean H column density (accounting for the 
presence of helium) is written as: 

(N H (t)) ~ 6.8 10 21 g(t)/ H atoms cm" 2 . 
where g(t) is the gas fraction. A galaxy with g ss 0.2 has 
< Nh >^ 1.410 21 atom cm~ 2 , in good agreement with 
the observational value for late-type galaxies (GRV87). 
/h is a parameter which allows us to increase the H col- 
umn density and is close to unity for normal spirals. The 
extinction curve here depends on the gas metallicity Z(t) 
according to power law interpolations based on the Solar 
Neighbourhood and the LMC and SMC with s = 1.35 for 
A < 2000 A and s = 1.6 for A > 2000 A. The derived ex- 
tinction curves depend strongly on the adopted geometry, 
here we restricted ourselves to the "slab" geometry; keep- 
ing in mind that it might introduce some uncertainties in 
the obscuration curves. It would be interesting to compare 
this with more detailed models of extinction curves (see 
for example Ferrara et al. 1999) but this is beyond the 
scope of this present work. 

Finally, the observed flux F* x is then related to the 
emitted (intrinsic) flux F£ by: 



(3) F t 



— = IQ-° AA ^ 



(5) 



The bolometric luminosity released at IR/submm wave- 
lengths is then given by: 



L IR (t) = J F A *(i)(l - 10-° AA ^)d\. 



(6) 
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2.3.2. Dust emission 



In order to compute the IR/submm emission spectra, 
one needs to redistribute the total IR luminosity Ljr over 
the whole range of wavelengths. This is achieved by using a 
dust model of the ISM. We assume the three-component 
(Polycyclic Aromatic Hydrocarbons, Very Small Grains 
and Big Grains) dust model developed by Desert et al. 
(1990) according to Maffei's method (Maffei 1994) which 
uses the observational correlations of the IRAS flux ratios 
with Ljr. 

Three components are considered in Desert et al. model 
spectra: 

— the Big grains (BG). The are made of silicates and 
graphite and have sizes between 10 nm and 0.1 /x m and 
are almost in thermal equilibrium. They are reasonably 
described by a modified blackbody € v B v (Tbg) with 
emissivity e„ = v m (where the index 1 < m < 2. 

— Polycyclic aromatic hydrocarbons (PAHs) . Because these 
molecules are small in size (less than 1 nm) they never 
reach thermal equilibrium when excited by the UV/vis 
radiation. Their temperature fluctuates and can exceed 
the equilibrium temperature which produces the 12 /im 
excess and the bands at 3.3, 6.2, 7.7 and 11.3 /im. 

— Very small grains (VSG) with sizes between 1 and 10 
nm also made of silicates and graphite. As PAHs they 
never reach thermal equilibrium and therefore their 
emission spectrum is much broader than the black body 
spectrum at a given temperature. 

The temperature Tbg of the Big Grains is given by 
the IRAS 60/100 colours provided we fix the value of the 
emissivity index m. The derived temperature is between 
15 and 50 K. Here we adopt the standard value of m = 1.5. 
Model spectra with m = 1.5 seems to cover all the range of 
observations in the submm (Guiderdoni et al. 1998). The 
contributions of each component is then calculated iter- 
atively from the 12/100, 25/100 and 60/100 ratios. The 
resulting spectra are computed from few /jm to several 
mm and evolve with the total IR luminosity L/^ such that 
galaxies with higher L/jj emit preferentialy at shorter A 
(Fig. |l]). Here we will adopt the model spectra of dust 
emission with total IR luminosity of Lm = 10 10 L Q which 
is typical of normal galaxies and corresponds to a temper- 
ature of BG of Tbg ~ 16-5 K. In spite of its limitations, 
this method is able to reproduce the FIR photometry of 
various individual galaxies (Maffei 1994). 



2.4-. Assumptions and inputs of the model 

We choose three classes of CSFR. The high SFR model 
(model H) choosen in order to be consistent with the 
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Fig. 1. Model spectra in the IR/submm for different total IR 
luminosities ranging from 10 6 to 10 14 Lq and an emissivity in- 
dex m = 1.5. 
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Fig. 2. The star formation rate as function of redshifts. The 
curves are the different star formation histories we input in 
our models: (solid line), (small dashes) and (long dashes) cor- 
responding to H, Li and L2 models. The data points are taken 
from M96 and F99 (filled squares): lower points and upper 
points respectively. Treyer et al. 1997 (empty triangle), Con- 
nolly et al. 1997 (empty hexagons). At high-z, the data are 
from Madau et al. 1998 (open triangles) and from S99 after 
correction for extinction (crosses). 
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recent estimations based on IR surveys (F99) and large 
ground-based survey of Lyman-Break galaxies (S99). The 
low SFR models L\ and Li are consistent with the UV/optical 
derived SFR with respectively a declining and a flat SFR 
at high redshifts (see table f ). Fig. || shows the three mod- 
els of CSFR we used, together with a compilation of the 
CSFR estimations derived from UV/vis and IR surveys. 
In our modelling of dust obscuration we used the HDF 
rest-frame attenuated 1500 A emissivity at z = 3 — 4 
(from S99 for models H and Li and from M96 for model 
L 2 ) and the local emissivity at 100 fim to determine the 
two adjustable parameters, Ty(z = 0) and a. 
In the case of the closed-box approximation 
0) the equations of chemical evolution (eq. 1) are solved 
assuming that at early times the stars form from some 
amount of gas whose initial gas density fraction is fixed 
by imposing that the present-day amount of HI gas is of 
the order of O s (0) «5x 10~ 4 (as estimated by Briggs & 
Rao 1993). In the case where the closed-box approxima- 
tion is relaxed, we will assume for simplicity a constant 
outflow tlg Ut — const. Note here that the ISM gas content 
of the galaxies is assumed to be mainly due to neutral hy- 
drogen gas HI; we neglect both the ionized and molecular 
forms of the gas (HII and H2). This approximation is rea- 
sonable at high redshifts, as measurements of gas content 
of Damped Lyman a systems show low ionization states 
(Viegas 1995) and low molecules abundances (Levshakov 
et al. 1992, Petitjean et al. 2000). Locally, this approxima- 
tion may not hold since the present-day late-type galaxies 
contain on average a non-negligible fraction of H 2 (Young 
& Scoville 1991) which would result in a factor of two un- 
certainty on Qg as or more if H2 contribution to the ISM 
gas is more important (Pfcnnigcr et al. 1994). 
We adopt a Salpeter IMF <fi(m) cx mT x with x = 1.35 for 
nid < m < m up and rrid = 0.1 M© and m up — 100 M . 
This IMF has been found to be more reliable in reproduc- 
ing the observed emissivities locally than the Scalo IMF 
(Lilly et al. 1996, M98). 

Finally, we will suppose to a first order that the results 
are independent of the cosmological parameters. Discus- 
sion of the effect of the cosmological parameters on the 
luminosity density can be found in (Casse et al. 1998). 



model 


< z < 1 


z > 2 


H 


consistent with F99 


consistent with S99 


Li 


consistent with M96 


consistent with S99 


U 


consistent with M96 


consistent with M96 



Table 1. Definition of the class of models describing the CSFR 
histories we input in our modelling. H model is chosen to be 
consistent with the mid-infrared CSFR inferred from observa- 
tions of ISO 15 /im selected galaxies in the CFRS at < z < 1 
by F99 and with the CSFR inferred from large ground-based 
survey of Lyman-break galaxies by S99 and corrected for ex- 
tinction at z > 2 . Models L\ and L2 both are consistent with 
the CSFR derived from the CFRS survey by M96 at < z < 1, 
at high-z model L\ is consistent with S99 while L2 is consistent 
with M98 




-2-10 1 2 3 

log X[ /x] 

Fig. 3. The rest-frame spectral energy distribution SED from 
UV/vis to FIR/submm at z = 0. The data points are from 
Ellis et al. (1996) for 0.44 /im ; Gardner et al. (1997) for 2.2 
pirn; and Soifer & Neugbauer (1991) for 25, 60 and 100 /im. 



3.1. The rest-frame comoving emissivity from 
UV/optical to IR/submm 



3. The results and comparison to observations 

Deep galaxy surveys have allowed the probing of rest- 
frame UV/optical to NIR emissivities of high and low red- 
shift galaxies. When combined with the CIB measurement 
it provides us with crucial constraints on the global his- 
tory of star formation. 



We are interested in the global evolution of the radi- 
ation from both the stellar and the dust component. This 
can be described by the comoving luminosity density p v 
as a function of redshift defined as the energy radiated per 
unit frequency per unit of co-moving volume at redshift z. 
The intrinsic stellar emissivity can be expressed in terms 
of the global star formation rate f2* 

pl = Pc f F^{t-t')hs'W (7) 

Jo 

where F£(t—t') is the stellar population spectrum as com- 
puted by our population synthesis code (see DGS for de- 
tails). 
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The spectral luminosity density at each redshift is com- 
puted by simply summing the observed (attenuated) stel- 
lar radiation and the dust radiation: p v (z) = pl x (z)+pf(z) 
where p* x is related to the intrinsic (unabsorbed) emissiv- 
ity by: p* oX = lO" 04 ^. 

We first fix the adjustable parameters of the models a 
and Ty, by fits to the observed luminosity density at f500 
A at z = 3 — 4 and the local to the local 100 ^m emis- 
sivity respectively. The best-fit solutions are obtained for 
the following values of the parameters (ry, a): (1.3, 1.4), 
(1.2, 1.4) and (1.2, 0.8). The corresponding values of the 
extinction factor Ay varies from Ay — 0.45 locally to 
Ay = 2 — 2.5 at z = 3 — 4 for model H. In the case 
of model L<x, Ay is almost constant, varying from 0.45 
locally to 0.53 at z = 3 — 4. Fig. || shows the resulting 
z = co-moving spectral luminosity density from UV/vis 
to IR/submm wavelengths, on which are superimposed the 
data. Our models with a Salpeter IMF give overall good fit 
to the present-day stellar light from UV/vis to IR. Except 
for the 2.2 pm emissivity which seems to be overestimated 
by if -type models, our models are able reproduce the data 
not only at 100 pm (which is expected since it is used as 
an input) but also at A = 0.44, 2.2, 25, 60 pm, whatever 
the chosen CSFR history. This result shows that local co- 
moving emissivities are not sensitive to the CSFR and 
therefore are not suitable to constrain the star formation 
rate history. We also compute the resulting spectral lumi- 
nosity density in the case where GRV87 prescription for 
T\(t) is used and found no difference with models which 
use our simple prescription (eq. 2); the two SEDs are al- 
most undistinguishablc. 

We use for the determination of dust emission a model 
with Lift = 10 10 -L©. To take into account the uncertain- 
ties in the FIR part of the SED, we run our models with 
a higher IR luminosity L]r > 1O 12 L which is typical of 
the Ultra-Luminous Infrared Galaxies (ULIRGs). As ex- 
pected, we found that the peak emission in Far-Infrared 
(100 pm) is shifted toward shorter wavelengths, leading to 
a poorer fit to the data. 



3.2. The Extragalactic Background Radiation 

The second observational constraint is set by the ex- 
tragalactic background light (EBL). Cosmological impli- 
cations of the CIB on the CSFR have been discussed by 
several groups, for example M98, Dwek et al. (1998), Gis- 
pert et al. (2000). Here we use our models to derive the 
intensity of the background. The intensity l v of EBL re- 
sults from integrating the emission coming from various 
extragalactic objects on the line of sight over the cosmic 
time. 

The intensity I v at a given time t is given by the following: 



We here ignore the contribution of AGNs in heating the 
dust. Predictions of our models for the intensity of the 
background radiation in the optical and IR/submm [] are 
shown in Fig. || The overall spectrum seems to be well re- 
produced by all the models regardless of their very differ- 
ent SFR histories. In particular, the shape of the submm 
part of the extragalactic background does not discriminate 
between the models, and is therefore not a very crucial 
probe of the star formation history. Instead, the CIB data 
point at 140 pm is very discriminant. Only mid-infrared 
derived CSFR at low-z (model H) is able to correctly 
match this point while UV-derived low- z CSFR (Li and 
L2 models) , falls considerably short of the COBE detec- 
tions at A = 140 pm. Moreover, the CIB is less sensitive 
to the CSFR shape at high-z. For example, had we cho- 
sen a model similar to //-model at low-z but with a de- 
clining SFR at high-z, we would obtain an accpetable fit 
to the 140 pm data. This result suggests that the bulk 
of the energy in the FIR/submm is due to moderately 
obscured normal galaxies lying at low and moderate red- 
shifts (0 < z < 1) and not by a population of heavily ex- 
tinguished star-forming regions located at very high red- 
shifts, as previously claimed (Blain et al. 1998, hereafter 
B98). We also show the predicted intensity of the back- 
ground radiation when the GRV87 prescription for the 
optical depth t\ is used (eq. 4). As we can see there is no 
strong difference between the two prescriptions. 
It is interesting to compare our result with other works. 
PFH have made use of many observational inputs (from 
QSOs absorption line surveys, optical imaging, redshift 
surveys, and FIRAS/DIRBE cosmic infrared background) 
to derive solutions for the cosmic histories of star forma- 
tion rate. They found a global star formation rate which is 
higher than the UV-derived one, in good agreement with 
our result. Although their three solutions do not fit the 
CIB at 140 pm, their solutions are actually significantly 
below this data point (see their fig. 11). Our result in 
the other hand, is in disagreement with previous works in 
which it is claimed that the best CSFR solution is higher 
and flatter at high redshifts than the optically derived 
form (B98, Gispert et al. 2000). Note that these authors 
used a different approach than ours and did not use the 
point at 140 pm to discriminate between the models, as 
in the present work. 



4. Cosmological history of metals 

In order to get a consistent picture of the evolution of 
the universe, one should relate the CSFR to the metal en- 
richment of the universe and its stellar content. To achieve 
this goal we will use the CSFR model that account for all 



c f°° 



(8) 



1 Recent tentative of CIB detection at 2.2 and 3.5 pm has 
been reported by Gorjian et al. (2000), our best-fit model is 
consistent with their estimated intensity at 2.2 and 3.5 pm 
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log \[fj,m] 



Fig. 4. The predicted diffuse backgrounds in the optical and 
FIR/submm as compared to current limits and detections. The 
coding of the curves are same as in fig 1. Solid hexagons re- 
sult from the faint galaxy counts at 3600 - 22000 A (Pozzetti 
et al. 1998). The empty square results from galaxy counts at 
2000 A (Armand et al. 1994). The Empty hexagon is the detec- 
tion at 3.5 ^m from COBE/DIRBE by Dwek & Arendt (1998). 
The filled squares correspond to the COBE/DIRBE residuals 
(Hauser 1996) only 140 (jm and 240 /im signals are firm de- 
tections (Hauser et al. 1998), all other points are upper limits. 
Empty triangles are the 100, 140 and 240 /im estimations from 
Lagache et al. (2000). The Filled triangles are the estimations 
of the CIB at 850 and 450 (jm from repsectively Blain et al. 
1999) and Smail et al. (1997) updated by Ivison et al. (1999). 
The Empty squares are the lower limit to the background at 
7 and 15 ^ m from Altieri et al. ( 1999). The empty losanges 
are the lower limit to the 15 and an upper limit to the 
0.912 fim background by Vogel et al. (1995). The zigzag is the 
COBE/FIRAS residuals detected by Puget et al. (1996) and re- 
visited by Guiderdoni et al. (1997). The lower panel shows the 
predicted diffuse background intensity obtained using GRV87 



available data from UV/vis to FIR/submm wavelengths 
to compute the enrichment history of the universe. 

4-1. Chemical enrichment of the ISM 

Our knowledge of chemical evolution was based on 
local observations (mainly the solar neighbourhood) until 
the outcome of high-z observations of damped Lya sys- 
tems, which allowed us to draw some conclusion about 
galaxy evolution in terms of the cosmic evolution of the 
mean metallicity and cool neutral gas. There are some in- 
dications that these DLAs represent the progenitors of 
present-day galaxies, although the nature of these ob- 
jects is still a matter of debate. Moreover, it is not yet 
clear whether the mean HI weighted metallicity derived 
from current observations of DLAs is representative or 
not of the global metal enrichment of the Universe. The 
picture that emerges from the DLAs observations is that 
i) the neutral gas content is decreasing from z ~ 3 to 
the present-epoch, corresponding to the conversion of gas 
into stars and that ii) the metallicity is decreasing with 
increasing redshift, with very low metallicities at high red- 
shifts, which would suggest that these systems are in their 
early phase of enrichment. 

Figure || shows the mean abundance of metals in the inter- 
stellar medium Zism a s function of redshifts compared to 
the metallicity at high redshift derived from the observed 
abundance of Zn in damped Lya systems (Pettini et al. 
1997), which should agree if DLA and field galaxies have 
the same metal enrichment. The uncertainty in the data is 
still large though, and no attempt has been made to cor- 
rect for the observational biases which could potentially 
alter our interpretation of the DLAs metal abundances 
(Prantzos & Boissier 2000 and reference therein) but see 
Savaglio 2000. We find that in the closed-box assumption, 
our best-fit to the CIB at 140 fxm model clearly fails to 
reproduce the data; it overpredicts the mean metallicity 
by a factor of^3at2<z<3. This is less than that 
found by B98 with their submm-derived CSFR (a factor 
as high as 5). Our best-fit model also produces a mean 
metallicity Z(0) which is ~ 2 times larger than the metal- 
licity in the local universe if we assume it to be solar. 
We have therefore relaxed the closed-box approximation 
and allowed for some outflowing of metal-enriched gas 
from the galaxies. As we can see from figure 0, the agree- 
ment with the data becomes better than in the closed-box 
models. If we further assume that the ejected metals are 
responsible for IGM pollution, we can use our model to 
estimate the amount of metals present in the IGM. 



4-2. Chemical enrichment of the IGM 

The detections of CIV and Si IV absorption lines in 
the Lya forest as reported (Tytler et al. 1995, Songaila 
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& Cowie 1996, Songaila (1997) can potentially give us in- 
formation about the early enrichment of the intergalactic 
medium (IGM) with metals (Gnedin & Ostriker 1997). If 
we assume that the IGM metal content can well be rep- 
resented by the metallicity inferred from Lyev clouds, we 
may compare our prediction of the density of metals to 
the observations of CIV and Si IV absorption lines in the 
Lya forest. Assuming that ejected metal-enriched gas is 
uniformely distributed in the intergalactic medium, we es- 
timate the mean co-moving density of metals in the IGM 
by integrating Clf ut over the time. We find that the den- 
sity of heavy-elements at z — is = 7.7 x 10~ 5 (i.e 
p m = 5.4 x 10 6 M Q Mpc- 3 ). At z w 3.2, Q^ M - 6 x 1(T 7 
in agreement within the error bar with the data point of 
fi£f M - 8 x 1CT 7 reported by Songaila (1997). This result 
supports the general belief that the chemical enrichment 
of the IGM may be due to the outflowing of metals from 
galaxies, although this scenario is still a matter of debate 
(see Gnedin 1997 for an alternative scenario). 
It has been pointed out that if no differences exist be- 
tween field and clusters galaxies, then clusters of galaxies 
can provide us with information about the chemical en- 
richment history of the universe (Renzini 1997). We there- 
fore estimate the local IGM mean metallicity Zjgm(0) = 
Vt^ M /Vl b and find it to be of the order of 1(T 3 . Using the 
low value of D/H ( r) W ~ 5.3) by Buries & Tytler (1998) 
(fit, ~ 0.077) and assuming a solar metallicty of 0.016, 
this value yields to a global mean metallicity of 6% solar 
which is ~ 4 times smaller than the metallicity derived 
from X-ray observations of the intra-cluster medium ICM 
(~ 0.25Z Q ). As suggested by Renzini (1997), this cluster- 
field difference may be due to either cluster-related pro- 
cesses such as stripping or ram pressure or to a flatter 
IMF in clusters relative to field galaxies. Note that if we 
use the high value of D/H ps 10~ 4 (e.g. Lemoine et al. 
1999 for a review) yielding fi^ ps 0.02, we would obtain a 
global metallicity of 0.2 solar which is close to the ICM 
metallicity. 

Regarding the large uncertainties in the observations and 
the limitations of our models (a constant metal-enriched 
gas ejection is certainly a crude approximation) , we should 
be cautious before drawing any conclusion regarding the 
global enrichment of the universe. 



4-3. Stellar content of the universe 

Our models can also provide the global evolution of 
the stellar content. Fig. || shows the inferred evolution of 
the stellar content as a function of the cosmic epoch for 
each of the models of CSFR . The estimate of the current 
2 = !1,(0) ~ 0.004 - 0.008 comes from Briggs (1997). 
We can see that the best-fit model to the 140 fim data 
point predicts a present-day stellar mass density larger 
by at least a factor of two than the current value. This is 



-3 




redshift 

Fig. 5. Evolution with redshift of the mean comoving ISM 
metallicity as predicted by the best-fit model (solid line). The 
dotted line curve corresponds to the case of closed box evolu- 
tion. The data points correspond to the metallicity as inferred 
from the DLAs Pettini et al. (1997). The point at z = 0.77 is 
from Boisse et al. (1998). 



very important as it shows that models able to produce 
enough energy in the FIR part of the CIB lead to an excess 
of stars as compared to the present-day estimates. Note 
that the FIR/submm derived CSFR of B98 produce values 
of f2*(0) at least a factor of two higher than that which 
our models produce; this is due to their higher CSFR at 
high redshifts. The current f2*(0) estimate is still uncer- 
tain though. There are several sources of errors that can 
bias the inferred f2*(Q), such as the uncertainties in the 
derived luminosity functions parameters $* and L*; on 
the estimation of the contribution of starlight from low 
surface brightness (LSB) galaxies; in the adopted stellar 
M/L ratio of disks, spheroids and irregulars. However all 
these effects are not very important and cannot account 
for the discrepancy. More exotic solutions should there- 
fore be invoked. Several solutions are possible, first it is 
possible that some stars (extragalactic stars, machos) are 
missed in the census of the stellar budget but it is unlikely 
that we missed that significantly. 

Second, contamination by AGN causes problems that are 
more important than we think. We know from IRAS and 
SCUBA observations and recently from X-ray Chandra 
observations (Barger et al. 2000) that a population of dust 
enshrouded AGNs do exist (Sanders and Mirabel 1996). 
However it is still unclear how important this contamina- 
tion is as it is difficult to discriminate between a pure star- 
burst and a dusty AGN in the object identification pro- 
cess. Therefore, it is highly probable that a large fraction 
of the CIB may be due to the dust enshrouded AGN pop- 



R. Sadat, B. Guiderdoni and J. Silk: Cosmological history of stars and metals 



9 



ulation. AGN contamination also complicates the deriva- 
tion of the star formation rate history from the observed 
luminosities even in the multi-wavelengths surveys (as the 
one derived by F99), although it is less critical than in the 
UV wavelengths. 

Third, our adopted Salpeter IMF is responsible for such a 
discrepancy. A IMF biased against low-mass stars, would 
reduce the number of long-lived stars and thus attenuate 
the discrepancy. Evidence of a shallow mass function be- 
low 1 M has been observed in the Galactic disk (Gould et 
al. 1996) and more recently in the Galactic bulge (Zoccali 
et al. 2000). A heavily biased IMF has also been suggested 
from microlensing data (Chabrier 1999). 




1 2 3 4 5 



redshift 

Fig. 6. Evolution as a function of redshift of of the co-moving 
stellar mass density. Coding of the models are as in fig. [| The 
horizontal lines define the data of the total mass density of the 
stars at the present epoch as estimated by Briggs (1997) 



5. Conclusions 

We have used a series of self-consistent chemo-spectro- 
photometric models which take into account the metallic- 
ity and dust obscuration effects to study the global evo- 
lution of spectrophotometric properties together with the 
stellar, chemical and radiation content of the galaxies av- 
eraged over a large co-moving volume, independently of 
their individual details and properties. 

• We find that our models with a Salpeter IMF repro- 
duce fairly well the observed present-day comoving lumi- 
nosity density at A = 0.44/im and 2.2 /urn although mid- 
infrared derived CSFR (if-model) appears to overproduce 



the light at A = 2.2^m. At wavelengths dominated by 
the dust emission (at A = 25, 60 and 100 /im), we find 
good agreement with the observations independent of the 
choosen CSFR. 

•We show that the point at 140 fim is crucial for con- 
straining the CSFR history and seems to favour a higher 
comoving star formation rate at low redshifts than the 
UV/optically derived one; however the shape at high red- 
shifts of the CSFR is not well constrained by the CIB 
measurements, contrary to previous claims. 

• We suggest that the main contribution to the bulk 
of the CIB and more specifically to the energy at 140 /jm 
is due to normal galaxies, lying at low and moderate red- 
shifts, and not to distant and dusty galaxies as previously 
suggested. 

• Concerning the history of the metal enrichment, we 
find that in the frame of closed-box evolution, our best-fit 
model to the CIB at 140 /im overpredicts the metallicity 
as observed in DLAs. A better agreement is obtained in 
the case of a model with some outflowing metal-enriched 
gas. 

• Assuming that the ejected metal-enriched material 
is responsible for the IGM metallicity we derive, the IGM 
metal content and found it to be consistent within the 
error bar with the metals in the Lya forest. However, 
we cannot derived conclusions regarding the IGM metal- 
enrichment given the large uncertainties in the data and 
the crudness of our metal enrichment treatment. 

• Our best-fit model to the CIB at A = 140 over- 
poduces the present-day stellar mass density f2*(0) by a 
factor of at least ~ 2 as compared to the current value. If 
the CIB at A = 140 /im is correct, then this discrepancy 
and the excess of light at A = 2.2 /Ltm could be explained 
by our choice of Salpeter IMF. This would then be an 
indication in favour of an IMF which is biased toward 
massive stars, opening a new debate on the universality of 
the initial mass function. However, one may keep in mind 
that our calculations have been done under the assump- 
tion that the FIR/submm background is due to a stellar 
component only; we neglect any contribution of AGNs in 
dust heating, which is probably not true. 
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